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ABSTRACT 

Model color magnitude diagrams of low-metallicity globular clusters usually show a deficit of 
hot evolved stars with respect to observations. We investigate quantitatively the impact of such 
modelling inaccuracies on the significance of star formation history reconstructions obtained from 
optical integrated spectra. To do so, we analyse the sample of spectra of galactic globular clusters 
of Schiavon et al. with STECKMAP (Ocvirk et al.) and the stellar population models Vazdekis 
et al. and Bruzual & Chariot, and focus on the reconstructed stellar age distributions. Firstly, 
we show that background/foreground contamination correlates with E(B-V), which allows us to 
define a clean subsample of uncontaminated GCs, on the basis of a E(B-V) filtering. 

We then identify a "confusion zone" where fake young bursts of star formation pop up in 
the star formation history although the observed population is genuinely old. These artifacts 
appear for 70-100% of cases depending on the population model used, and contribute up to 12% 
of the light in the optical. Their correlation with the horizontal branch ratio indicates that 
the confusion is driven by HB morphology: red horizontal branch clusters are well fitted by old 
stellar population models while those with a blue HB require an additional hot component. The 
confusion zone extends over [Fe/H]= [—2,-1.2], although we lack the data to probe extreme 
high and low metallicity regimes. As a consequence, any young starburst superimposed on 
an old stellar population in this metallicity range could be regarded as a modeling artifact, if 
it weighs less than 12% of the optical light, and if no emission lines typical of an HII region 
are present. This work also provides a practical method for constraining horizontal branch 
morphology from high signal to noise integrated light spectroscopy in the optical. This will allow 
post-AGB evolution studies in a range of environments and at distances where resolving stellar 
populations is impossible with current and planned telescopes. 

Subject headings: galaxies: stellar content, galaxies: formation, galaxies: evolution, techniques: spec- 
trosopic, stars: horizontal branch, globular clusters: general 



1. Introduction 

Spectroscopy is a major tool for the study of 
galaxy formation and evolution. The analysis of 
the spectrum of galaxies provides constraints on 
their star formation activity, past and present. 
However the accuracy of these constraints depends 
on the quality of the data, but also critically on the 
quality of the models used to interpret the s e spec - 
tra, as shown for example in Koleva et al.l ( 2008 ). 
Since spectral models of stellar populations rely on 



stellar evolution models, it is vital the latter mod- 
els be able to account for most of the observed 
stellar evolution stages, or at least the most lu- 
minous ones (e.g. turnoff, red giant branch and 
horizontal branch stars). While thermally puls- 
ing asymptotic giant branch (TP-AGB hereafter) 
stars are fundamental for NIR spectral modelling, 
hot evolved stars are crucial in the UV for interme- 
diate and old populations. The latter have been 
proposed as t he origin of the U V- upturn of ellip- 
tical galaxies (|Han et al.ll2007t Iy3 12008). In this 
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paper, we will focus on these hot horizontal branch 
(hereafter HB) stars. They are fairly luminous for 
their mass: while weighing only about 0.5 Mq, 
they have a typical absolute magnitude My=-0.7 
((Coxll2000h . Their He-burning core tends to be 
hotter with decreasing metallicity, maki ng the HB 



of low -Z galactic globular clusters bluer (|Lee et al 



1994|), although the morphology of the HB also de- 



pends on a second parameter, so far elusive. 

Accurate modelling of the HB morphology 
of GCs is a difficult problem. For instance, it 
i s strongly s ensitive to mass-loss on the RGB 
( Potter 2008). Moreover, accurate fitting of blue 
HB in GCs usually requires a fine tuning of He 
abundanc e dispe rsion, as shown in NGC2808 by 



Lee et al.l (|2005l ). This clashes with the gen- 
eral treatment of spectral synthesis for galaxies, 
where a single stellar population (hereafter SSP) is 
chemically homogeneous. In effect, the parameter 
space allowed by popular full spectrum SSP mod- 
els (as opposed to those pre dict ing indices) such 
Bruzual & Charlotl (|2003l) or iLe Borgne et al ' 



a.s 



(|2004l ) is strictly (age, metallicity), and the abun- 
dances are solar-scaled, with no intrinsical disper- 
sion. With only these parameters, it is expected 
that the isochrones adopted in these SSP models 
are not appropriate for the whole range of HB 
morph ology. Indeed, Fig. 7 of Bruzual k, Charlotl 
(|2003l ) shows that the 'Padova 1994' and 'Padova 
2000' tracks do not fit well the color-magnitude 
diagrams (hereafter CMD) of GCs with blue HB. 
In particular the HB of the models does not ex- 
tend enough towards the blue. These very tracks 
are nonetheless used for the computation of spec- 
tral SSP models. What can be the impact of 
this lack of hot evolved stars when interpret- 
ing the spectra of stellar populations ? Their 
effect can be readily seen in the in tegrated col- 
ors of GCs (|Smith &: Straderl I2007L Moreover, 
recent works indicate that space observations in 
the GALEX NUV and FUV are required to dif- 
ferent iate between blue HB stars and young ho t 
stars (jKavirai et al.ll2007HSchawinski et al.ll2007j ). 



However, although GALEX has been successful 
at securing data for a fair fraction of the sky, a 
large corpus of past, present and ongoing extra- 
gala ctic studies still r elie o n optical spectra onl y 
/e.g. iReichardt et al l (l200ll):IPanter et al.[ (120031) 
Heavens et all (120041) ; rToieiro et alT 



It is now well known that HB morphology 
makes age-dating stellar populations from op- 
tical spectra only troublesome , usin g indices 
(|Beaslev et al.l l2000i: IPuzia et al.l 120051). but also 
using full spectrum fitting techn iques (|Koleva et al 
2008; ISharina &: Davoustl 12009): the age inferred 
spectroscopically can be several Gyr younger than 
that measured from the turnoff. The effect on star 
formation history (hereafter SFH) reconstructions 
could be even more dramatic: do the missing hot 
stars sho w up as an additiona l fake burst as sug- 
gested in ICoelho et al.l (2009), or do they simply 
shift the whole age distribution to younger ages ? 
At the moment, no studies have been published 
specifically on the effect of HB morphology on 
the significance of SFH reconstructions from full 
spectrum fitting in the optical. 

In this paper, we propose to investigate quan- 
titatively this very issue. To do so, we first need 
a sample of objects with well-known, simple star 
formation histories, and a variety of well con- 
strained HB morphologies. From this point of 
view, GCs are advantageous since they are mostly 
believed to arise from a single star formation 
episode, and display a whole range of HB mor- 
phologies. We will focus on the library of inte- 



grated spectra of galactic GCs of iSchiavon et al . 
(|20n5h with STECKMAP (jOcvirk et al.ll2006bl lal). 



a Bayesian method for reconstructing the star 
formation history of galaxies using full spec- 
trum fitting, and a collection of popular sin- 
gle stell ar population (SSP) models based on 
MILES dSanchez-Blazauez et al.l 120061 ) and also 
Bruzual & Charlotl (|2003E For each GC, the 
CMD is available from the literature, so that we 
can relate HB morphology to possible departures 
of the SFH reconstructions from a unique star 
formation burst. We recall that the aim of the 
paper is not to propose a new spectral model for 
the sample of GCs: it is well established that 
interpreting HB morphology in globular clusters 
requires more parameters tha n the simple SSP 
population models used here dBedin et al. 2004; 



Koleva et al.l |2009h : Michielsen et al.l (|200 



Ferraro et al ]|2004tlNorris1l2004l ). Instead, our goal 
is to understand how the unknown HB morphol- 
ogy in galaxies may complicate the interpretation 
of their spectra. That's why our setup reproduces 
the prevailing regime in which most extragalactic 
stellar populations studies in optical integrated 
light spectroscopy are carried out and therefore 
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uses SSPs or combinations of SSPs. The plan 
of the paper is as follows: Sec. [2] presents the 
data, the SSP models and the method, and Sec. 
presents the recovered stellar age distributions 
and their dependence on HB morphology. In Sec. 
IH wc discuss the implications of our results for 
spectroscopic extragalactic studies. 

2. Methodology 
2.1. Data 

We consider the sample of integrated spec- 
tra of galactic glob ular clusters published by 
Schiavon et al. ( 2005 ^1. It consists of 40 Galac- 
tic globular clusters, obtained with the Blanco 
4 m telescope and the R-C spectrograph at the 
Cerro Tololo Inter- American Observatory. The 
spectra cover the range ~ 3350 — 6430 A with 
« 3.1 A (FWHM) resolution. The spectroscopic 
observations and data reduction were designed to 
integrate the full projected area within the cluster 
core radii in order to properly sample the light 
from stars in all relevant evolutionary stages. The 
signal to noise values of the flux-calibrated spectra 
range from 50 to 240 per A at 4000 A and from 
125 to 500 per A at 5000 A. The selected targets 
span a wide range of cluster parameters, including 
metallicity, horizontal-branch morphology, Galac- 
tic coordinates, Galactocentric distance, and con- 
centration. We dropped NGC6569, NGC6388 and 
NGC6441, which have rather peculiar HB mor- 
phology. We are then left with 37 GCs. The pa- 
rameters relevant to this study ar e listed in Tab 
1, and are taken from Table 1 of Schiavon et al 



(|2005l ). In order to characterize HB morphology, 
we use the Horizon tal Branch Ratio (HBR) indi- 
cator, defined as in Lee et al. ( 19941 ): 



HBR = (B - R)/(B + V + R) , 



(1) 



where B, V, R are the number of blue, vari- 
able and red HB stars, respectively. GCs with 
a blue HB have HBR « 1, while GCs with a 
red HB have HBR w -1. For NGC5946 and 
NGC628 4, we assign HBR = 1 (no value is pro- 



I 2005f n. 



Since the GCs we 
stu dy here have masses of the ord er 10 5 — 10 6 M Q 
(see lChernoff fc Diorgovs ki (1989) for instance for 



mass determinations), and the scanned area en- 
closes a large fraction of the mass, we consider that 
our results should not be affected by stochastic 
fluctuations related to the fin iteness of their stellar 
popul ations. According to lLancon k, Mouhcine 
(l2000h . the flux fluctuations for our sample should 
be smaller than a few percents. 

2.2. Models 

We consider two popular SSP models. The 
first one (in preparation, hereafter Vazdekis 
et all) is based on the MILES stellar li- 



brary (jSanchez-B lazqucz ct al 



isochrones from Girardi et al 



120061) and the 
( 2000f) . The other 



one is lBruzual fc Charlotl (|2003h (he reafter BC03) 
based on the STELIB stellar library (|Le Borgne et al 
l2003h. and using the "Padova 1994" tracks 



(|Bertelli et al.ll994h . The corresponding isochrones 
for [Fe/H]= [-0.7,-1.3,-1.7] are shown in Fig. 
[TJ Note that throughout the paper we will use the 
notation [Fe/H]= log ln (Z /Z ffi ), and Z© = 0.0189 
(| Anders fc Grevessdll989h . 

The isochrones are very similar up to the begin- 
ning of the RGB. In the latter phase, the Girardi 
2000 set is bluer, especially for [Fe/H]=-1.3. Then 
both sets develop a horizontal branch which ex- 
tends further to the blue with increasing age and 
decreasing metallicity. 

The spectral models we use thus include the HB 
contribution by construction. However, although 
the predicted position of the BHB in the CMD is 
rather satisfactory in these isochrone sets, the pre- 
dicted number of stars in this phase is very small. 
For instance, a synthetic CMD obtai ned with IAC- 
star0 in the default setup using the lGirardi et al 



(|2000t ) isochrones for a [Fe/H]=-1.7, 17 Gyr old 
population, produces a horizontal branch ratio 
HBR = 0.2, while most GCs of our sample have 
HBR = 1 at the same metallicity and an age of 
about 10 Gyr. In effect, we checked carefully the 
spectral evolution of the low-Z models from 10 to 
20 Gyr in the optical and they show very little 
sign of the increasing BHB contribution with age: 
for the Vazdekis et al. models, a very slight blu- 
ening and deepening of the Balmer lines (about 
1%) can be seen at these extreme ages, however 
for the BC03 models we see nothing at all. In any 



1 http: / / www.noao.edu/ggclib / dcscription.html 



2 http: / / www.iac.es / galeria / vazdekis / vazdekis_models_ssp_seds .html 
3 http: / /iac-star. iac.es/iac-star/ 
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Fig. 1. — Color-magnitude diagra ms of the isochrones used for the Vazdekis et al. models (jGirardi et al 
(|2000l ). black solid line) and BC03 (jBertelh et all (|l994l ). grey dashed line) for 3 metallicities [Fe/H] = [-0.7,- 
1.3,-1.7]. For each panel, the sequence spans the [4-20] Gyr age range by logarithmic increments of 0.1 
dex. 



case the contribution of the BHB to the optical 
spectrum even in 20 Gyr models is much smaller 
than the difference in GCs spectra induced by dif- 
ferences in HB morphology for clusters of simi- 
lar rne_tall_idty J _as shown for instance by Fig. 3 
of ISchiavon et al.1 (|2005h . So even if the BHB is 
present in the isochrones, it is legitimate to say 
that there is almost no BHB contribution to the 
spectra in the optical within our set of SSP models. 

Both models cover equally well the wavelength 
range of the data. We consider the metallicity 
set [Fe/H] = [-1.7,-0.7,-0.4,0,0.4]. Although 
BC03 covers a wider age range, we restrict our- 
selves to SSP sequences of age=[10 8 — 2.10 10 ] yr 
for both models, divided in 40 age bins. We in- 
clude these extreme ages even though they are 
larger th an the WMAP5 age of the universe of 
13.7 Gyr (H inshaw et al.ll2009l ). because, as shown 
by Fig. [TJ really blue horizontal branches de- 
velop only at ages larger than 12 Gyr in both 
sets of isochrones. We did not use PEGASE-HR 
(|Le Borgne et al.ll2004h coupled w ith the ELODIE 
3.1 library (jPrugniel et alJ 120071 ). although it of- 
fers the highest spectral resolution (R= 10000), be- 
cause it does not cover the A < 3900 A domain. 

2.3. Interpretation software: STECKMAP 

In order to analyse the observed spectra with 
SSP models as a reference, we use the interpre- 



tation code STECKMAP (|Ocvirk et al.ll2006bl fal). 
The SSP models are used in a flux-normalized 
form (i.e. the average flux of ever y SSP is 1 ), 



rather than mass-normaliz e d, as i n lOcvirk et al 
(|2006bh and lOcvirk et all (|2006ah . In this set- 
ting, STECKMAP returns the most likely stel- 
lar age distribution (hereafter SAD) correspond- 
ing to the analysed spectrum, i.e. the luminous 
weight of each model component as a function 
of age, thereby providing an analog of the star 
formation historjff], giving insight into the lumi- 
nous balance between various stellar components 
seen in the spectrum. One important aspect of 
STECKMAP is that it is a regularized method, 
i.e. the various unknowns (here SAD, age-Z re- 
lation, and line-of-sight velocity distribution) are 
forced to be smooth in order to avoid exaggerat- 
edly oscillating solutions. The smoothing kernel 
is laplacian for the SAD ( i.e. L = D2 using the 
notations of Ocvirk et al.l ( 2006af )). and we chose 
/ix = 1, which gives robust results. 

Besides, we require the age-Z relation to be fiat 
(i.e. mono-metallic population), and implement 
this constraint by using a finite difference opera- 
tor L = Di — diag 2 [— 1,1] and a large smooth- 
ing parameter = 10 5 . The smoothing kernel 



4 The star formation history can be obtained by multiplying 
the stellar age distribution by the mass-to-light ratio of 
each bin divided by the time spanned by each time bin 
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Fig. 2.— Spectrum of NGC1851 (black) and best fitting model using STECKMAP and the Vazdekis et al. 
SSPs (red). The green line represents the non-parametric transmission curve as explained in Sec. I2.3I Bad 
pixels and some residuals of sky subtraction have been masked, as well the very blue and the red end of the 
spectrum beyond 6200 A. 



for the line-of-sight velocity distribution is lapla- 
cian as for the SAD, and the smoothing parame- 
ter was set to \i = 10 2 , which produces a smooth 
gaussian LOSVD for the whole sample. In order 
to deal with possible flux calibration errors, we 
mutliply the model by a smooth non-parametric 
transmission curve, representing the instrumen- 
tal response multiplied by the interstellar extinc- 
tion. This curve has 30 nodes spread uniformally 
along the wavelength range, and the transmission 
curve is obtained by spline-interpolating between 
the nodes. The latter are treated as additional 
parameters and adjusted during the minimization 
procedure. This continuum matching technique is 
similar in essence t o the multiplicative polyn omial 
used by NBursts (jChilingarian et al.|[2007h . and 
will be detailed in a forthcoming paper. 

Finally, the first guess for the solution consists 
of a flat, constant SAD. However, in some occa- 
sions, we used a family of SADs with a unique 
gaussian burst centered on a random age to probe 
the existence of local minima (see Sec. I3.2I) . 

3. Results 

We now present the results obtained by apply- 
ing STECKMAP with the described setup to the 
sample of GCs. A very good fit to the data is ob- 
tained in all cases, as shown for NGC1851 in Fig. 



2. 



3.1. 2 families of stellar age distributions 

We find that the recovered SADs fall in two 
classes, as illustrated by Fig. [3l 

• type I: genuinely old SAD, with a single 



bump consisting only of > 
panel). 



10 



9.5 



yr (left 



• type II: two-component SAD, with most of 
the light provided by an old population, and 
a small light fraction (up to 12%) coming 
from the youngest bins (right panel). 

While the first family of SADs is what one 
would actually expect for GCs, the second one 
comes as a surprise, since GCs are usually ex- 
pected to consist of stars formed in one single 
short star formation event. We already pointed 
out that hot evolved stars such as BHB stars, if 
unaccounted for in the SSP models, could show up 
in the spectra by producing abnormally deep blue 
Balmer lines. In such cases, the model SSPs will 
be redder and display weaker Balmer lines than 
the observed spectrum. As a consequence, our in- 
terpretation code STECKMAP indicates that the 
best fitting combination of the SSP models pro- 
vided consists of an old SSP plus a very young 
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Fig. 3. — Stellar age distributions of NGC1851 and NGC6218 using the Vazdekis et al. models. 



SSP that mimics the HB stars unaccounted for by 
the model isochroncs. 

We now test the relevance of this interpreta- 
tion of type II SADs by investigating quantita- 
tively their connection with HB morphology. 

3.2. Hot fraction and horizontal branch 
morphology 

We define the hot fraction fhot as the sum of 
the luminous weights of the young (age < 10 85 
years) components. We recall here that the SAD 
is normalized (the sum of all weights is 1). In 
Fig. 2] we plot the hot fraction fhot versus the 
HBR for the whole sample. About half of the 
clusters have fhot = 0, which was our initial ex- 
pectation. Moreover, there seems to be a general 
trend of increasing fhot with HBR. About 10 of 
the GCs have HBR « 1 and fhot > 0, which sup- 
ports our interpretation for the type II SADs, as 
explained in Sec. 13.11 There are however out- 
liers with respect to this simple picture: 4 GCs 
with red HB have fhot > 0, and another 3-4 have 
a blue HB but fhot ~ 0. We checked that these 
were not secondary minima by using a family of 
single bursts populations centered on random ages 
as first guess. They all led to the same solution. 
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Fig. 5. — Hot fraction as a function of color excess, using the Vazdekis et al. models. Left: red HB clusters, 
right: blue HB clusters. 



3.2.1. Background/ foreground contamination 

The two objects in the top left part of Fig. Q] are 
NGC6553 ad NGC6528, which are located in the 
direction of the galactic center. Because of their 
location, they are definitely the most likely to be 
contaminated. 

Since we are interested in the effect of HB mor- 
phology on the spectra rather than contamina- 
tion, we wish to work on a sample free from fore- 
ground/background contamination. One could de- 
fine a subsample by removing the GCs too close to 
the galactic center direction. However, this does 
not provide a good filtering, since the likelihood of 
a contamination depends on the integrated stellar 
density inside the scanned volume and the expo- 
sure time, and hence also on the distance to the 
GC. We thus investigated the use of alternative 
indicators of contamination. 

Red HB morphology is usually well reproduced 
by the tracks used for stellar populati o ns sy n- 
thesis (see Fig. 7 of lBruzual fe Chariot! (120031 )): 
unlike for blue HB GCs, old SSP models suffice 
to reproduce the spectra of red HB GCs. Be- 
cause of this superior agreement with models, we 
can use them as calibrators: any deviation from 



type I SAD in a red HB GC should be caused 
by background/foreground contamination rather 
than inaccurate evolutionary models. Therefore, 
for these objects, fhot should trace contamina- 
tion. One could object that the large fhot found 
could also be a consequence of a strong contribu- 
tion from blue stra gglers (hereafter BS). Indeed, 



Zoccali et al.1 (|200ll ) shows a clear "blue plume" in 
the CMD of NGC6553. However, their kinemati- 
cal cleaning shows that only half of the stars of the 
plume are indeed cluster members and therefore 
genuine BS. The other half of the plume consists 
of young disk stars in the foreground, and is hotter 
and more luminous than the BSs. Hence, even if 
BSs certainly partially contribute to the blue light 
excess, the latter is dominated by foreground disk 
contamination. The high fhot found for NGC6553 
thus effectively traces contamination rather than 
abnormal BSs frequency. 

The left panel of Fig. [5] shows fhot as a function 
of color excess for GCs with HBR < 0. Indeed, 
most of the red HB GCs are well fitted by gen- 
uinely old SSPs, without the need for any addi- 
tional hot component. However, for E(B-V)> 0.4, 
fhot increases sharply, reaching a maximum for the 



7 



2 galactic center GCs NGC6553 and NGC6528. 
This correlation between fhot and extinction lends 
support to our interpretation. In effect, a correla- 
tion between background/foreground contamina- 
tion and extinction should naturally arise from the 
two following points: 

• if dust and stellar density correlates, then 
E(B-V) increases with the integrated stellar 
density along the line of sight to the GC. 

• young hot contaminating stars are more 
likely to live in dusty environments. 

It seems therefore that the color excess E(B- 
V) is a good proxy for the likelihood of fore- 
ground/background contamination, and the left 
panel of Fig. [5]tells us that GCs with E(B-V)< 0.4 
should be relatively free from contamination. In 
order to clean our sample from contaminated spec- 
tra, we will drop the high-extinction GCs, and 
keep only those with E(B-V)< 0.4. The right 
panel of Fig. [5] shows that, within this new defini- 
tion of the sample, there is no correlation between 
E(B-V) and f hot for the blue HB GCs either, and 
the points are randomly distributed. It also shows 
that the E(B-V) filtering only removes 4 BHB 
GCs, and leaves us with a large enough sample to 
carry on with our investigation. 



3.2.2. The f, 



hot 



HBR correlation 



The left panel of Fig. [6]shows the results for the 
cleaned sample. The correlation between fhot and 
blue HB morphology is now much clearer: spuri- 
ous bursts of recent star formation are very likely 
to appear in GCs with HBR > 0, while fi 10 t = 
for all the GCs with HBR < 0. The intensity of 
the spurious young burst can be as high as 10% of 
the total optical light of the cluster. The spurious 
recent bursts correlate with HB morphology. 

However, we find that 4 (out of 13) GCs with 
BHB have an almost fhot ~ 0. To check this last 
puzzling point and more generally the robustness 
of the fhot-HBR correlation we investigated the in- 
fluence of the theoretical error in the SSP models. 

3.3. Model dependence 

We repeated the analysis of the sample using 
the BC03 models. The results are shown in the 
right panel of Fig. [6l First of all, the correla- 
tion between fhot and HBR is recovered, and is 
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Table 1: Properties of the sample of GCs and re- 
sults. (1): NGC name, (2): color excess, (3): 
Metalliciy [Fe/H], (4): horizontal branch ratio, (5) 
Fraction of light in the hot component in %, com- 
puted from the SADs obtained with the Vazdekis 
models, (6): same as (5) with BC03 . Col umns 
(l)-(4) are taken from lSchiavon et aT (2005), and 
for (5)-(6), see Sec. El 
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Fig. 6. — Hot fraction as a function of horizontal branch morphology for the clean sample, using Vazdekis 
ct al. (left), and BC03 models (right). 



therefore a robust result of our study. However, 
the frequency of spurious young bursts is higher in 
general with this model, both in GCs with BHB 
and RHB. In particular, the low-Z GCs for which 
fhot = with the Vazdekis models, now all have a 
measurable fhot > 0. On the other hand, 2 of the 
red HB GCs, which previously had fhot = now 
have a non-zero hot fraction. 

The overall consistency of the results obtained 
using the 2 models indicates that, for the regime 
studied here (sub-solar metallicity, age > 6 Gyr), 
the systematic error introduced by the uncertain- 
ties in the modelling of the HB are larger than 
the other error sources (differences in the under- 
lying stellar libraries, accuracy of the stellar pa- 
rameters and coverage of the latter, interpolation 
procedures). 

3.4. Method dependence 



In lKoleva et al.l (|2008l ). the authors performed a 
similar experime nt, however with a differen t, para- 
metric method ([Chilingarian et al.l I2007I ). where 
the young SSP was replaced with a hot star spec- 
trum. They also found a correlation between the 
HBR and the luminous weight of the hot star, al- 
though much weaker than the one presented in 



this paper in Fig. [5] In a different type of exper- 
iment, NBURSTS and STECKMAP were found 
to give consistent results for luminosity-weighted 
ages and metallicities. Moreover, ULYSS (Koleva 
et al., in prep.), another non-parametric code, was 
also found to give similar results to STECKMAP 
in the reconstruction of the s tar formation histor y 
of dwarf elliptical galaxies (jKoleva et al. 2009). 
Because of this agreement, we suggest that the 
clearer f w-HBR correl a tion w e find here with 
respect to Koleva et al. ( 2008l ) is essentially the 
consequence of the wider wavelength range of 
the models we used. Indeed, the Vazdekis et 
al. models allow us to exploit the full cover- 
age of t he data into the b lue, down to 3500 
A. while iKoleva et al.l (l2008h used PEGASE-HR 
(|Le Borgne et al.ll2004l h which starts at 3900 A. 
This gives us good confidence in the fact that our 
results are robust to a change of method, as well as 
a change of SSP models, provided that the wave- 
length range are kept similar. 
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4. Discussion 

4.1. The confusion zone 

Wc identified a range in metallicity ([Fe/H] < 
— 1.2) where the analysis of the integrated spectra 
of stellar populations in the optical is highly likely 
(70—100% probability depending on the model) to 
lead to the detection of a spurious young burst of 
star formation, contributing up to 12% of the light. 
The luminosity of this spurious burst correlates 
with HB morphology, and increases towards bluer 
HB. This result is independent of the SSP model 
and of the method used. 

The consequence for studies of stellar popula- 
tions of galaxies is that any detection of a young 
burst of star formation superimposed on an old 
metal-poor population ([Fe/H] < —1.2) from anal- 
ysis of the optical spectrum is possibly spurious 
if the contribution of the young populations is 
smaller than 12% and no emission lines typical of 
HII regions are seen. Moreover, since galactic stel- 
lar populations are composite in general, smaller 
young fractions can still be the signature of BHB 
stars even if the old component has a luminosity- 
weighted metallicity larger than —1.2. For exam- 
ple, following our results, an old composite popu- 
lation, with 50% of mass in solar metallicity stars 
and 50% of mass in [Fe/H] = —1.2 stars could 
still display 5% of the light in a fake young burst. 
What exactly does it mean for extragalactic stud- 
ies ? In the following we discuss the implications 
of our result for spectroscopic studies of various 
objects in integrated light. 

4.2. Field and cluster HB stars 

So far we have considered exclusively GCs. It 
is not obvious that our result also applies to in- 
tegrated spectroscopy of field stars, which is the 
method of choice for extragalactic studies. 

Indeed, some stages of stellar evolution seem 
to differ depending on membership to a GC 
or to the field, suc h as metal-rich RR - Lyrae 



stars (Lavden et al. 19991 : Catelan et al. 120061: 



Matsunagal2007 ) . However, these differences seem 



to disappear at low Z (i.e. the regime of this work). 
Moreover, differences in the properties between 
field and GCs HB stars (other than RR Lyrae) 
have not been detected. It is not clear why some 
stages are affected by GC membership while oth- 



ers are not. Two processes at least are expected to 
be significantly more important in cluster environ- 
ments than in the field. These are binary stellar 
evolution (higher binarity in clusters), and bi- 
nary disruption (decr easing binary frequency), as 
sho wn for example in iGuhathakurta et al. (1998) 
and ISollima et al.l (|2007l ). As a matter of fact, 
blue HB morphology is in deed obs e rved in the 
field as well: the CMDs of iDolphinl (|2002l ) show 
that the dwarf galaxies Leo I and especially Ursa 
Minor have a well defined BHB. In the absence 
of strong contrary evidence, we consider that the 
HB of the GCs analysed here are representative of 
the HB morphologies of at least a fraction of the 
field stellar populations of galaxies. 

4.3. Metal-rich populations 

The highest metallicity object of our clean 
sample is NGC104, with [Fe/H]=-0.7. We were 
unable to use the higher metallicity part of the 
original sample because it happens to be also 
the most extincted, and thus the most contam- 
inated. However, upon visual inspection of the 
CMDs of NGC6528 and NGC6553 ([Fe/H]=-0.1 
and [Fe/H]=-0.2), no evidence for a blue HB is 
found, and these populations should therefore be 
well behaved with respect to star formation his- 
tory reconstructions. This extends our results to 
quasi-solar metallicities. At even higher metal- 
licities, we recall that spectroscopy of elliptical 
galaxies in the optical suggests sm all but recent 
star formation (Trager et al.ll2000h . in the form 
of a "frosting" of young stars, reminiscent of the 
artifacts we have observed for GCs with BHB. 
This is a very important aspect since the halt 
of star formation in massive galaxies, as inferred 
from their red colours, motivates a large amount 
of studies in the search for quenching mecha- 
nisms, either through hot shocks and the end of 
cold s t reams ( Birnboim &: Deke l 2003; K eres et al 



2005: Ocvirk et al. 12008). or AGN feedback 



([Hopkins et alj l2007t |Pe Lucia fc Blaizotl |2007[) . 
We note however that evidence for recent star for- 
matio n in early-type galaxie s is also found in IR 
data ( Lvubenova et all l2008h . which is little sen- 
sitive to any EHB/BHB contribution. To probe 
such high metallicities, an integrated spectrum 
of NGC6791, an old (8-9 Gyr), metal-rich (twice 
solar) open cluster, would be very useful. Inter- 
estingly, several extreme horizontal branch stars 
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(EHB), even hotter than usu al BHB stars, have 
been found in this obj ect (jKaluznv &; Udalski 
Il992t iLiebert et J] 1 19941 ) . although without any 
sizeable BHB counterpart. However, it is unlikely 
that this small population of very hot evolved 
stars can produce such strong artifacts in SFH 
reconstructions as we have seen from optical light. 
Their contribution peaks in the FUV and affects 
the optical light much less than the BHB studied 
here. As an alternative to NGC6791, one could 
study old metal-rich extragalactic GCs, although 
we would then miss the accurate photometric de- 
termination of their HB morphology from a CMD. 

4.4. M32 

The rightmost panel of Fig. 6 of Coelho et al. 
2009 shows that except for the nucleus, a large 
fraction (about « 50 %) of the light of M32 is pro- 
duced by a [Fe/H] < — 1.2 old stellar population. 
Our study shows that in this regime, a spurious 
young burst accounting for sa 5% is indeed likely 
to appear. However, in the nucleus, the underly- 
ing stellar population is metal-rich. Even though 
we do not constrain well the extent of the "confu- 
sion zone" in the high Z regime, the young compo- 
nent is stronger than the fake bursts we found: it 
accounts for 20 % of the light (in a similar wave- 
length range) while we found maximum spurious 
contributions of 10-12%. Hence, our result sup- 
ports the claims of Coelho et al. 2009, mainly 
that there has been a recent burst of star forma- 
tion in the center of M32, while the young compo- 
nent seen in the outskirts is also compatible with 
being an artifact due to the blue HB morphology 
of the system. 

5. Conclusions 

The analysis of int e grated spectra of galac- 
tic GCs ISchiavon et all (|2005l ) with STECKMAP 
yields two families of stellar age distributions: 



• Type I: genuinely old populations, with no 
contribution from components younger than 
yr. 
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Type II: old populations with a superim- 
posed hot component younger than 10 85 yr 
contributing a fraction fhot to the light of the 
GC. 



GCs with red HB and a type II SAD are subject 
to background/foreground contamination. Since 
the latter correlates with color excess E(B-V), we 
define a clean sample (with respect to contamina- 
tion) by removing GCs with E(B-V)> 0.4. In this 
new sample, red HB GCs all have type I SADs, 
while blue HB GCs have mostly type II SADs, 
with up to 12% of the light in the hot component. 
This is the signature of inaccuracies in the mod- 
elling of the horizontal branch: in the model tracks 
used in Vazdekis et al. and Bruzual &: Chariot! 
(J2003J), the blue HB is often too short compared to 
observations. As a consequence, the observed inte- 
grated spectrum contains an additional hot com- 
ponent with respect to the corresponding model 
SSP, which gives rise to the fhot-HBR correlation. 
This is the main result of the paper, and it is found 
with both SSP models Vazdekis et al. and BC03. 

We have thus identified a "confusion zone" , 
where spurious starbursts are very likely to appear 
in star formation history reconstructions, even- 
though the observed population is genuinely old. 
These fake starbursts are young (< 10 8 5 yr) and 
occur in 70-100% of cases depending on the SSP 
model, and can weight up to 12% of the light in 
the optical. This confusion is driven by HB mor- 
phology and is thus mainly a low-metallicity effect 
([Fe/H] < —1.2). In this regime, any minor young 
burst of star formation can be an artifact, es- 
pecially if no emission lines associated with star 
formation are seen. We expect this prediction to 
be independent of the interpretation method and 
models, as long as the evolutionary tracks used 
for the spectral synthesis are not tuned specifi- 
cally to account for BHB morphology. The ex- 
act weight allowed for the fake starburst depends 
on the wavelength range considered. It should in- 
crease for domains bluer than in this work, and 
decrease towards the red. 

The highest mctallicity of our clean sample 
is [Fe/H] =-0.7, which prevents us from properly 
probing the high-Z regime. We plan to do so by 
studying old metal-rich extragalactic GCs or an 
integrated spectrum of NGC6791. 

We wish to encourage the stellar evolution and 
extragalactic communities to consolidate the re- 
cent progresses made on HB morphology model- 
ing by including the relevant tracks into future 
spectral synthesis models such as Vazdekis et al.. 
However we note that accurate fitting of blue HB 
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morphology in GCs by fine tuning of He a bun- 
dance dispersion such as in iLee et all (|2005l ) in- 
troduces additional parameters, and possibly new 
degeneracies in spectrum interpretation. More- 
over, proper use of these tracks requires a stellar 
library, empirical or theoretical, for a range of non- 
solar abundance ratios, reliable everywhere on the 
CMD. Such models are stil l in their infancy, d e- 



spite improvements such as ICoelho et al.l (|2007fl 



This study is also an interesting sanity check 
for STECKMAP, which is found to behave very 
well. We propose that the exercise of analysing 
the E(B-V)-filtcrcd sample and checking the re- 
sults dependence on HBR as we have done be a 
test for all the star formation his tory reconstruc- 



tion methods, such as MOPED ([Reichardt et al 
20011). VESPA dToieiro et al.ll2007l). STARLIGHT 



(jCid Fernandes et al.ll2004) or lWild et alj (|2007l ). 
and urge other authors to perform the same ex- 
periment. 

Finally, the fact that STECKMAP associates 
BHB morphology with a type II SAD with a high 
success rate, shows that HB morphology can be 
determined in integrated light through full spec- 
trum fittin g, and in the optical domain, as was 
realised bv lSchiavon et al.l (|2004 ). This technique 
could thus become a valuable tool for extragalac- 
tic cluster studies: by removing the need for UV 
(i.e. space) observations and resolved populations 
(i.e. relatively nearby), it gives access to samples 
and environments otherwise unreachable. 
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